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Abstract. Multiple, sequential mergers are unavoidable in the hierarchical build-up picture of galaxies, in particular 
for the minor mergers that are frequent and highly likely to have occured several times for most present-day 
galaxies. However, the effect of repeated minor mergers on galactic structure and evolution has not been studied 
systematically so far. We present a numerical study of multiple, subsequent, minor galaxy mergers, with various 
mass ratios ranging from 4:1 to 50:1. The N-body simulations include gas dynamics and star formation. We study 
the morphological and kinematical properties of the remnants, and show that several so-called "minor" mergers 
can lead to the formation of elliptical-like galaxies that have global morphological and kinematical properties 
similar to that observed in real elliptical galaxies. The properties of these systems are compared with those of 
elliptical galaxies produced by the standard scenario of one single major merger. We thus show that repeated 
minor mergers can theoretically form elliptical galaxies without major mergers, and can be more frequent than 
major mergers, in particular at moderate redshift. This process must then have formed some elliptical galaxies 
seen today, and could in particular explain the high boxiness of massive ellipticals, and some fundamental relations 
observed in ellipticals. In addition, because repeated minor mergers, even at high mass ratios, destroy disks into 
spheroids, these results indicate that spiral galaxies cannot have grown only by a succession of minor mergers. 
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1. Introduction 

Major mergers between spiral galaxies of similar masses 
are known to form elliptical-like galaxies. The remnants 
of such violent events present an r^^^ r adial mass pro- 
file, as observed in eUiptical galaxies ( de Vaucouleurd 
19771) , and are mainly pressure-su pported. This is the case 



for equal-ma ss mergers (e.g. iBarnes fc Hernguistl 11991 



Barneslll992[ ). but more generally for m ass ratios ranging 
from 1 : 1 to 3:1 and even 4:1 , studied b y 'Bendo fc B arnesI 
( 2000l ) , ICretton et all (|200ll ) , iNaab fc Burkcr t (2003) and 
in Bournaud, Jog & Combes (2005b, hereafter paper I). 
The merger remnant tends to be boxy for 1:1 mergers, 
while 3:1 and 4:1 mergers mostly result in disky elliptical 
galaxies - that are not disk galaxies, but elliptical galaxies 
with a disky isophotal shape. An outer disk-like compo- 
nent appears in particu lar for mass ratios of 3:1 and higher 
( Naab fc Truiillol I2OO6I ). It can be present even in 1:1 or 
2:1 merger remnants but is far from being the main contri- 



bution to the total stellar mass. The disk component be- 
comes more massive with increasing mass ratios. Beyond 
4:1, mergers do not form elliptical galaxies any longer, but 
early-type disk-dominated systems (Bournaud, Combes & 
Jog 2004, and paper I). They are usually called "minor" 
mergers, but we have suggested in these earlier works that 
one should distinguish between "intermediate" mergers 
(between 4:1 and around 10:1) that form SO- like galaxie^, 
and the real "minor" mergers with ratios larger than 10:1, 
where the remnants can be classified as (disturbed) spiral 
galaxies (Quinn, Hernquist & FuUagar 1993; Velazquez & 
White 1999; Walker, Mihos & Hernquist 1996). 

Existing studies of multiple mergers have focused 
on nearly-simultane ous mergers of several galaxies of 
comp arable masses (|Wen fc Hernguistl Il994 il996.: Bekkil 
200lh . This typically corresponds to the collapse of a 
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^ These systems have a disk-like profile, but a massive bulge 
and their support is dominated by pressure instead of rotation , 
unlike spirals, m aking them rather SO-like l|jog fc Chitre|[2002l : 
iBournaud et al.ii2004 ') 
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compact group into one sing le giant elliptical galaxy. 
Weil fc HernauistI (|l994 Il996li have shown that the el- 
liptical galaxies produced by these events have some kine- 
matical properties that differ from the remnants of pair 
mergers, a nd provide a better explanation for observed 
properties. iBekkil (j200ll ) has studied the starbursts occur- 
ring when such events involve gas-rich late-type spirals, 
and the possible connection with ULIRGs. There is obser- 
vational evidence that some gala xies are remnants of the 
collapse of compact groups (e.g., Borne et al. 2Qq3) or at 
leas t of the simultaneous me rging of more than 2 galax- 
ies (|Taniguchi k, Shiovalll998l ). However, these studies are 
limited to situations where the different galaxies merge 
nearly at the same time, and have comparable masses, so 
that the merger of only two of them would already have 
resulted in the formation of an elliptical galaxy. 

In this paper, we study a different process: the mul- 
tiple, sequential mergers of the intermediate and minor 
types - hereafter both called "minor" for simplicity - i.e. 
the mass ratios that do not form elliptical-like galaxies af- 
ter one single merging event. This corresponds to a given 
spiral galaxy that gradually grows by merging with smaller 
companions. The mergers are not assumed to occur at the 
same time, but one after the other, over a total timescale 
of a few Gyrs. We show that a few successive mergers 
with mass ratios between 5:1 and 10:1 lead to the grad- 
ual transformation of spiral galaxies into elliptical-like sys- 
tems. Under the effects of the sequential mergers, spiral 
galaxies become earlier-type spirals, then lenticular SO- 
like systems, and finally spheroidal objects with global 
morphological and kinematical properties similar to ob- 
served elliptical galaxies. This new scenario for the for- 
mation of elliptical galaxies is compared with the result 
of the "standard" binary major merger scenario, although 
the differences are not striking. We also show that even 
the growth of a spiral merging with 50:1 dwarf compan- 
ions also leads to the progressive destruction of the spiral 
disk and the evolution towards SO- like and more frequently 
elliptical-like galaxies, if the initial galaxy is assumed to 
significantly increase its mass by such mergers with dwarf 
companions. Repeated minor mergers occurring alone can 
then in theory be an efficient process to form a large num- 
ber of elliptical galaxies. However, we discuss that in re- 
ality other processes can reduce this efficiency, in partic- 
ular to maintain massive disks in spiral galaxies. At the 
same time, some elliptical galaxies can likely have formed 
via multiple minor mergers, in particular the boxy ones 
and/or those formed at low redshift when minor mergers 
largely dominate major ones. 

Section 2 contains the details of N-body simulations 
and analysis of the results. In Sect. 3, we analyze the prop- 
erties of the multiple merger remnants as a function of 
the mass ratios and the number of mergers, and study the 
differences with or similarities to the elliptical galaxies re- 
sulting from major mergers. The implications for galaxy 
evolution are studied in Sect. 4, and Sect. 5 contains a 
brief summary of the results of this paper. 



2. N-body simulations of multiple galaxy mergers 

2.1. Code description 

We have used the ^N-body FFT code of 

Bournaud fc Combed ( 2003f l. The gravitational field 
is computed with a resolution of 400pc. We used 4 x 10^ 
particles for the most massive galaxy. The number of par- 
ticles used for the other galaxy is proportional to its mass. 
Star formation and time-de pendent stellar mass-loss ar e 
implemented as described in lBournaud fc CombesI (|2002h . 

The star formation r ate is coinpute d according to the 
generalized Schmidt law (ISchmidtjll959h : the local star for- 
mation rate is assumed to be proportional to j-ig^ , where 

is the local two-dimensional density of gas. We chose 
= 1.5, as sug gested b y the observational results of 
KennicuttI \im §) (see also iBoissier et al.ll2006f ). The dis- 
sipative dynamics of the ISM has been accounted for by 
the sticky-particles s c heme described in Appendix A of 
Bournaud fc Combed (|2002t ). In this paper we employ elas- 
ticity parameters /3t=/3r-=0.7. 



2.2. Physical model for colliding galaxies 

Each galaxy is initially made up of a stellar and gaseous 
disk, a spherical bulge and a spherical dark halo. The vis- 
ible mass of the main (target) galaxy is 2 x 10^^ Mq. Its 



stellar disk is a Toomre ( 19641 ) disk of radial scalelength 
5 kpc, truncated at 15 kpc. Gas is distributed in a Toomre 
disk of scale-length 15 kpc and radius 30 kpc. The bulge 
and dark halos are Plummer spheres of radial scalelengths 
3 kpc and 30 kpc respectively. The dark halo is truncated 
at 70 kpc. The bulge-to-total mass ratio is 0.17 (bulge- 
to-disk: 0.2), and the dark-to-visible mass ratio inside the 
stellar disk radius is 0.7. This corresponds to a stellar disk 
of 14 X 10^° Mq, gas disk of 2.5 x 10^° M0, a bulge of 
3.3 X lO^" Mq, and a dark halo of 11.7 x 10" M© (up 
to r = 70 kpc). The number of particles is 1.3 x 10^ for 
stars, 0.93 x 10^ for gas, and 1.76 x 10^ for dark matter. 

The initial velo cities of particles for each co mponents 
are computed as in Bournaud k. Combed ( 2003f l. The ini- 
tial value of the Toomre parameter is Q = 1.7 over the 
whole disk. The gas mass fraction in the disk is 15%. In 
this paper, we wish to compare successive minor mergers 
to binary major mergers formed under the same condi- 
tions, so we do not vary the gas fraction, except for the 
50:1 mergers presented later in Sect. 14.1. "21 

The companion galaxies have a radial distribution of 
matter scaled by the square root of their stellar mass. For 
instance a 4:1 companion has a disk that is half the size, 
as is the bulge, etc.. Gas mass fraction, bulge-to-disk and 
halo-to-disk mass ratios are kept constant to avoid a too 
large number of free parameters. The free parameters are 
then: 

— the number of merging companions 

— the mass ratio of each companion 

— the orbital parameters of each companion 
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2.3. Orbital parameters 

The orbital parameters of a merging pair of galaxies have 
been described e.g. by IDuc et aL (I2OOOI) . These authors 
defined two angles 9 and $ for each galactic disk. 9 is 
the inclination of the orbital plane vs the primary disk 
plane = for a prograde encounter in the disk plane, 
9 = 180 for a retrograde encounter in the disk plane). $ 
is the azimuthal position of the secondary disk spin axis. 
We have run series of simulations with $ fixed to 30 de- 
grees and a value of 9 of either 30 or 180+30 degrees, 
for each companion (i.e. inclinations of 30 degrees but or- 
bits either prograde or retrograde). The first encounter is 
prograde with a corotating companion, the second one is 
retrograde, with a counter-rotating companion, etc.. We 
fixed the encounter velocity to 150 km and impact 
parameter to 45 kpc, both computed at an infinite dis- 
tance, neglecting the action of dynamical friction before 
the beginning of the simulatior|^. Since the velocity at an 
infinite distance is larger than zero, the initial orbits are 
not bound but hyperbolic with a total energy larger than 
zero, and dynamical friction dissipates the relative kinetic 
of the galaxies, eventually leading to a merger. 

The choice of an inclination of 30 degrees is made be- 
cause it is near the statistically average value of the incli- 
nation angle in spherical geometry. Assuming an isotropic 
distribution of the initial orbital planes of galaxies, the 
probability of an inclination i is p{i) oc cos(i), and the 
average inclination < i >— J ip{i) di ~ 30deg. This 
average choice is then representative of orbits that are 
neither coplanar nor polar. The inclination has not been 
varied in this paper, except for sequences of 10:1 merg- 
ers as described in Sect. 14.1.11 It is known from paper I 
and other studies of major mergers that the properties of 
a binary merger remnant depend more on the mass ratio 
than the inclination and other orbital parameters; the test 
in Sect. 14.1.11 confirms that the result of a minor merger 
sequence depends on the total merged mass more than 
the orbits, even if changing the orbit does result in some 
variations of the final properties. 

In the reference case of one binary equal-mass merger, 
choosing one orbit orientation (prograde or retrograde) 
may not be representative; we thus simulated the two cases 
and kept the average properties in the following study: we 
still call this "run 1x1 : 1" for simplicity, and the properties 
of a merger remnants depend more on the mass ratio itself 
then the orbital parameters (e.g., paper I and other ref- 
erences in the Introduction). Similarly, for the 3:1 cases, 
we ran two sequences (prograde-retograde-prograde and 
retrograde-prograde-retrograde) and kept the average re- 
sults. For higher mass ratios, we discuss larger number of 
mergers so that starting with the first orbit as prograde 
or retrograde has less influence - the influence of the or- 
bit orientation is anyway smaller than the influence of 



the mass ratio itself (paper I for binary mergers, see also 
Sect. I4.1TT] in this paper for multiple mergers). 

2.4. Temporal sequences of mergers 

The nomenclature for each run is as follows: NxM: 1 indi- 
cates that the galaxy has undergone N mergers with com- 
panions of mass 1/M of the target galaxy mass. 
We have simulated mass ratios M of 1, 3, 5, 7, 10, 15, and 
N varying from 1 to M (by steps of 2 for M=15). 

In sequences with mass ratios up to 10:1, the mergers 
are separated in time by At = 800 Myr. This means that 
each interloper will reach its pericenter at an instant com- 
puted to be At after the pericenter of the previous ones. 
In studies of binary mergers over the same range of mass 
ratios (paper I) we found that this timescale ensures that 
each merger is relatively relaxe^l before the following one 
occurs - even if the merging/ relaxation occurs somewhat 
more rapidly for 1:1 cases than 10:1 cases. A time interval 
of 2 Gyrs or more would ensure a more complete relax- 
ation of the each merger before the next one occurs, but 
this would lead to total durations longer than the Hubble 
time. Actually, multiple mergers are rather expected to 
occur in groups or dense cosmological structures, with a 
total duration of the merger sequence of a few Gyrs and 
a somewhat short time interval between mergers. 

As for 15:1 mass ratios, the injection of fifteen com- 
panions each 800 Myr would represent a total duration of 
12 Gyr, which is highly unrealistic, in particular with the 
adopted mass and gas fraction of the initial target galaxy. 
We then chose to reduce the time interval between merg- 
ers to At = 400 Myr for this 15:1 mass ratio. The first 
merger is only partially relaxed when the second occurs, 
but the mergers are still subsequent, i.e. the nuclei of the 
first companion is merged when the second companion 
reaches the outer gas disk at r = 30 kpc. 

The time interval separating mergers is not varied in 
the present paper. We do not find strong differences be- 
tween the 10:1 {At = 800) and 15:1 {At = 400) merger 
sequences in the following, suggesting that it has only a 
minor influence on the main properties of the remnant 
galaxies, at least as long as the mergers are sequential 
and not simultaneous. Detailed properties like the orbital 
structure could however vary with the time mtcrval sepa- 
rating the merging events (e.g. lAthanassoulal2005h : in par- 
ticular, in the present paper we study subsequent mergers 
and do not compare to simultaneous mergers. 

In each run, the final system is evolved 1.5 Gyr after 
the last merger, to ensure full relaxation before any anal- 
ysis is performed. This way, the 5x10:1 simulation cor- 
responds to five mergers of mass ratio 10:1, separated in 
time by 800 Myr between their respective pericenter pas- 
sages, followed by a 1.5 Gyr isolated evolution to ensure 
the relaxation of the the final remnant before its prop- 



^ These values should be compared to the stellar disk radius 
of 15 kpc and circular velocity 180 km at the edge of this 
disk, to be applied to smaller/larger galaxies 



In particular, this 800 Myr interval after the first pericenter 
ensures that the nuclei from the first merger have coalesced 
when the second interaction is at its pericenter. 
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erties are analysed. In the the 10x10: 1 simulation, there 
are 9 intervals of 800 Myr between the mergers and a final 
1.5 Gyr-long relaxation of the remnant. 

We also ran a simulation with 50 dwarf compan- 
ions of mass ratio 50:1, and initial gas mass fraction of 
30%. This simulation, where the mergers cannot really be 
subsequent, but overlap over time, is described later in 
Sect. 14.1.21 We do not vary the orbital parameters beyond 
alternating prograde and retrograde orbits in the main 
part of the paper, but we will describe and discuss the in- 
fluence of orbital parameters in 10x10 : 1 merger sequences 
in Sect. HXn 

2.5. Analysis of the results 

Each remnant of a (sequence of) mergers is analyzed, 
once relaxed, to derive several physical quantities. We 
build projected column-density maps on which we define a 
"25th- isophote", its semi-major-axis being the R25 radius 
in the following, commonly called the "optical radius" in 
observations. The conversion of mass density into luminos- 
ity is calibrated so that the initial target spiral galaxy, seen 
face-on, has a central surface magnitude ^=21.7 for its 
disk component (not counting the additional l uminosity 
from t he bulge) . This is to be consistent with the FreemanI 
( 1970[ ) relation, where this central magnitude is primarly 
observed in B-band, but we do not explicitly assume any 
wavelength in the analysis of our simulations. We assume 
that the mass-luminosity conversion factor does not vary 
over the duration of our simulations, because we start with 
15% of gas so that the new stars will represent at most 
15% of the masflfl- Most stars are older than the merging 
events, hence have a reasonably constant mass-to-light ra- 
tio. The 25th isophote is found to include, on average, 82% 
of the stellar mass of the merger remants in our simula- 
tions. Once this reference radius R25 is defined, we can 
measure several fundamental parameters of the multiple 
merger remants: 

— the flattening E. We first derive the "face-on" projec- 
tion, defined as that minimizing the apparent flatten- 
ing of the 25th isophote. We then construct projected 
maps for ten "edge-on" projections, and measure the 
ellipticity e(r) = 10 x (1 — h/a) where the axis ratio 
6/a is obtained from an ellipse-fitting model. We finally 
keep the flattening E as the average value of e(r) over 
the [0.55i?25; -R25] range for the different projections. 
This definition is similar to that used in Bournaud 
et al. 2005b, and roug hly equivalent to that used by 
Naab fc BurkertI (|2003l) . 



the Sersic index n is determined as the best fit over the 
same [0.55i?25; ^^25] radial range, averaged for three or- 
thogonal projections (not restricted to "edge-on" pro- 
jections as was the case when deriving the flattening 
E). 



* actually most new stars formed in mergers lie in the central 
kpc, not around -R25, which results in an even lower contami- 
nation at this radius 



— the kinematical parameter V/ a. When building pro- 
jected images we also compute the velocity distribution 
along each line-of-sight, and then derive the average 
velocity V and velocity dispersion a for each one. We 
then define the V/a parameter as the column-density- 
weighted average (equivalent to luminosity-weighted 
assuming a constant mass-to-luminosity ratio) over the 
[0.55i?25; ^25] radial range, and the result is averaged 
over three orthogonal projections. The V/a parameter 
is not measured along a slit, but on all the pixels of a 
first momentum map within the chosen radial range, 
each pixel being weighted by its intensity in the asso- 
ciated density map. The choice to exclude the central 
regions (r < 0.55i?25) from this kinematical analysis 
was made to ensure that a disk galaxy will have a 
V/a parameter significantly larger than 1. Indeed this 
lower bound is always larger than 5 kpc in our merger 
remnants, which is larger than the bulges of spirals and 
SOs (even the large and massive bulges of intermediate- 
mass merger remnants in Bournaud et al. 2004). This 
way the presence of a massive rotating component will 
result in a F/cr > 1, while a value olV/a < 1 can- 
not relate to the bulge of a spiral galaxy, and indicates 
elliptical-like kinematics at large scales. 

— the boxiness 04. When measuring E from various pro- 
jections, we also compute the boxiness 04 (r) as a func- 
tion of the radius, and keep the average value 04 over 
the range [0.55i?25; ^25]- The choice of 0.55i?25 as 
the lower bound ensures that, when a massive disk 
is present, we really measure its boxiness or diskiness, 
not that of the bulge (see paper I). 

3. Formation of elliptical galaxies in successive 
mergers 

3.1. First detailed example 

We show in Fig. [T] the morphology and kinematics of the 
1x10: 1 to 4x10: 1 merger remnants. After one merger of 
this mass ratio, the system resembles an SO galaxy. It is 
made up of a massive disk, thickened by the interaction, 
and has large velocity dispersions but still with V/a > 1. 
After three or four 10:1 mergers, the remnant does not re- 
semble a disk galaxy any longer: it is a spheroidal galaxy, 
with maximal flattening E5-E6 (but less in most pro- 
jections). The isophotes have only a low residual diski- 
ness, and the Sersic index n ~ 3 shows that the radial 
luminosity profile is now closer to an "Jji/4" profile (de 
Vaucouleurs 1964) than exponential. The morphology of 
this system is hence typical of an elliptical galaxy. At the 
same time, the kinematics is dominated by velocity disper- 
sions with V/a = 0.6 after four mergers, similar to what is 
found in remnants of major 3:1-2:1 mergers (paper I and 
references in the introduction). 

The remnant of the succession of four 10:1 mergers 
has thus both the morphology and the kinematics of an 
elliptical-like galaxy, while each of these mergers would 
individually result in an early-type disk galaxy. After the 
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4x10 : 1 merger sequence, the system has the properties of 
a disky eUiptical galaxy where a significant residual rota- 
tion is still observed. With an increased number of minor 
mergers it can become a boxier, slower-rotating elliptical 
galaxy (see. Fig [5|). Also note in Fig. 1 that the resid- 
ual rotation axis also tends to become misaligned to the 
apparent flattening direction with increasing number of 
mergers. 

3.2. Morphology and kinematics of the multiple 
merger remnants 

The sequences of mergers with various mass ratios have 
been systematically analyzed. The fundamental parame- 
ters E, n, V/a, as defined in Sect. 2, are plotted in Figs.[2l 
[3] and m for each sequence. The results are given as a func- 
tion of the number of mergers converted into the "merged 
mass": for instance, when the initial galaxy has merged 
with 3 companions, each of them having a 5:1 mass ratio, 
this so-called merged mass is 1.6 (1 for the main initial 
galaxy and 0.2 for each companion). Using this definition 
allows us to directly compare sequences with various mass 
ratios. We recall that each multiple merger remnant is an- 
alyzed after 1.5 Gyr of evolution to make sure that the sys- 
tem is reasonably relaxed. We do not analyze the poorly 
relaxed systems between the occurrences of two successive 
mergers or the on-going merger properties in this paper. 

The first noticeable result is that once the merged mass 
is 2.0 (i.e., the total mass of the companions equal the 
mass of the initial galaxy, as in a 1:1 major merger), the 
system has the properties of a slow-rotating elliptical-like 
galaxy, be the companions massive (like 1:1 and 3:1) or 
much more minor (like 15:1). Indeed, all the final systems 
after these sequences of mergers have: 

— a spheroidal shape with an axis ratio larger than 0.5 
(class E4 or at most E5) when observed "edge-on". 
The shape of the systems is even rounder with other 
lines-of-sight. This is measured in the outer parts of the 
galaxies (up to R25), so that the systems cannot simply 
consist of spheroidal bulges surrounded by thin disks, 
but are really spheroid-dominated galaxies - if a thin 
outer disk is present, it lies beyond the 25th isophote 
and/or does not dominate the mass distribution. 

— a Sersic index of 4-5 on average, thus being much 
closer to the empirical law (de Vaucouleurs 
1964) than to the exponential profile of a disk galaxy. 
Here again the measurement is not restricted to the 
central regions but relates to the radial profile up to 
the 25th isophote. 

— a kinematical parameter V/a smaller than 0.3 (on av- 
erage 0.2), and are hence not rotationally-supported. 
Random motion dominates in these multiple merger 
remnants as in real elliptical galaxies, and are even 
larger than in the remnant of one single 3:1 merger 
(which is usually already considered to be an elliptical- 
like system). 
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Fig. 2. Evolution of the apparent flattening of the relaxed 
merger remnants seen edge-on, along sequences of merg- 
ers with various mass ratios. All systems evolve from the 
initial thin disk to spheroidal galaxies. The flattening ap- 
pears to depend on the merged mass, more than the mass 
ratio of the mergers: for instance 6 mergers of mass ratio 
10:1 or 3 mergers of mass ratio 5:1 (that both correspond 
to a merged mass of 1.6) both lead to the formation of an 
E5 system. Only the 15:1 systematically forms flatter sys- 
tems, but the difference remains small and several mergers 
of this type still form spheroidal galaxies. The values for 
1:1 and 3:1 mass ratios are averaged between prograde 
and retrograde orbits as detailed in the text. The "orbital 
scatter" shows the (Icr) dispersion of the result for the 
10x10: 1 sequence when orbital parameters are varied (as 
detailed in Sect. l4TT|) . 

There are differences from one case to the other, but 
they are not larger than the variations from one line-of- 
sight to another one for each given system, and much be- 
low the differences between standard elliptical and spiral 
galaxies. The main observable properties, in projection, of 
these remnants from multiple mergers are typical of ellip- 
tical galaxies, even if each merger is really minor (like 10:1 
or 15:1) and would by itself alone have kept the system in 
the form of a (disturbed) spiral galaxy. 

A merged mass of 2.0 is not required to form an 
elliptical-like galaxy. This is known to be the case for the 
"major" unequal-mass 2:1 and 3:1 binary mergers. A set of 
criteria to consider that a merger remnant is an elliptical 
galaxy can be: 

— a flattening of E6 or less when observed edge-on, i.e. 
thicker than any disk galaxy. 

— a kinetic energy mainly in the form of random motion, 
i.e. V/a < 1. 

— a Sersic index larger than 3, i.e. the luminosity proflle 
is closer to elliptical-like than to an exponential disk - 
an this cannot relate only to a central bulge since the 
inner regions are not included in our measurement. 

Thus, elliptical-like systems are obtained as soon as the 
merged mass is 1.3-1.4, i.e. the total mass of the compan- 






Fig. 1. Comparison of the initial spiral galaxy (here seen after 1 billion year of isolated evolution) and the relaxed 
merger remnants of one to four 10:1 mergers. The density maps (left) are in logscale. The velocity and dispersion fields 
are shown in linear scale, with color code ranging from -150 (black) to +150 km (white) for the velocity maps, and 
from (black) to 200 km s^^ (white) for the velocity dispersion. The contour intervals are 35 km s^^ in the velocity 
and dispersion fields. The system is viewed under the projection that maximizes its apparent flattening parameter E 
at each timestep; other projections make the system appear rounder, with generally lower rotational velocites. 



ions is at least 30% of that of the initial galaxy. This is 
in agreement with binary mergers up to 3:4 or 4:1 already 
producing elliptical galaxies (Bournaud et al. 2005b ~ po- 
tentially with faint outer disks: Naab & Trujillo 2006), 
but this also indicates than two mergers with mass ratio 
5:1 or three mergers with mass ratio 10:1 are enough to 
transform a spiral galaxy into an elliptical galaxy. These 



systems then resemble the low-mass elliptical galaxies that 
have a higher degree of residual rotation than the giant el- 
liptical galaxies (e.g., Naab & Burkert 2003). An increas- 
ing number of mergers finally keeps on evolving them into 
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Fig. 3. Evolution of the V/cr parameter along sequences 
of mergers of various mass ratios. All systems evolve from 
the initial rotating disk to remnants supported by increas- 
ing velocity dispersion. The kinematics of these multiple 
merger remnants appears to depend on the merged mass 
more than the mass ratio of the mergers. The velocity dis- 
persion becomes larger than the rotation velocity when 
the merged-mass is larger than ~1.3, after for instance 
one single 3:1 merger or three mergers of mass ratio 10:1. 
Very slow-rotating systems are obtained once the mass of 
the system has doubled via the succession of mergers, as 
after one single equal-mass merger. The "orbital scatter" 
is as in Fig. [2] (see also Sect. 14.1. 1|) . 
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Fig. 4. Evolution of the Sersic index n in the 
[0.55i?25; radial range. Starting from n ^ 1 for the 
initial disk, its value increasing while the luminosity profile 
becomes more typical of an elliptical galaxy. The "orbital 
scatter" is as in Fig. [2] (see also Sect. I4.1.T|) . 

more concentrated and roundeJl systems, together with a 
still decreasing degree of rotation. 

Thus, multiple minor mergers can form elliptical galax- 
ies that look like both the real observed ellipticals and the 



^ these systems are rounder regarding their global flattening 
£, but this does not imply that the diski-/boxiness is washed 
out (see Sect. 3.3). 



major merger remnants, from their shape (flattening), ra- 
dial profile (Sersic index) and kinematics iyja). These 
elliptical-like galaxies are never flatter than E7, like the 
real ellipticals and the major merger remnants. This is 
because the increase of the equatorial velocity dispersion 
is not achieved without a large increase of the vertical dis- 
persions, too, which makes the final system much thicker 
than its progenitors. 

With the parameters studied so far, no major system- 
atic difference exists compared to major mergers, provided 
that the total merged mass is the same. This can be in- 
terpreted as follows: the energy dissipated by dynamical 
friction is to first order equal to the relative kinetic energy 
of the merging galaxies (i.e. their mass and relative veloc- 
ities). The merging of an equal-mass pair thus releases the 
same energy through dynamical friction as ten subsequent 
10:1 mergers. The kinetic energy dissipated by dynami- 
cal friction heats the stellar systems. The same amount of 
heating is thus obtained in the two cases (binary 1:1 vs. se- 
quence of ten 10:1 merger), which explains the similar V jo 
final values. Kinematical heating in galaxy mergers is gen- 
erally close to isotropic, thus the different cases also lead 
to the same amount of vertical heating of the initial disk, 
which explains the similar flattening of the merger rem- 
nants of various origins. Finally, the rather similar Sersic 
indexes suggest that the final degree of relaxation is sim- 
ilar, even though it is reached through gradual stages. A 
detailed study of the dynamics and stellar orbits during 
on-going mergers will however be needed to fully under- 
stand the similarities and differences between the violent 
relaxation of major mergers and the step-by-step relax- 
ation through multiple minor mergers. 

3.3. Boxiness of elliptical galaxies 

Elliptical galaxies are often classified according to var- 
ious parameters. A fundamental one, that relates to 
their formation history, is the diskiness/boxiness param- 
eter 04, which quantifies the 4th Fourier componeni|f| in 
the a zimuthal decomposition of t heir projected isophotes 



fe.g-. lBender k Moellenhofilll987l ). Two classes of ellipti 



cal galaxies are generally distinguished: 

— the "disky" elliptical galaxies (04 > 0) are not disk 
galaxies but have an isophotal shape recalling that of 
disks, although at a much lower degree. These are gen- 
erally low mass elliptical galaxies, with a significant 
amount of rotation (V/a up to 0.5-1). In binary merger 
scenarios, these are mainly produced by 3:1-4:1 merg- 
ers. 

— the "boxy" elliptical galaxies (04 < 0), whose isophotes 
have a distorsion that recalls a box (or a rectan- 
gle) compared to perfect ellipses . These are generally 
higher mass elliptical galaxies, with less rotation {V/a 
generally limited to 0.1 or 0.2). In binary merger sce- 
narios, they result mainly from equal-mass 1:1 mergers 



^ We deflne 04 as the coefficient of the cos (4S) component 
in the azimuthal decomposition of the isophotes. 
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Fig. 5. Evolution of the 04 parameter versus the merged 
mass along several merger sequences. The elliptical-like 
galaxies formed with merged masses up to ~ 1.6 have a 
disky isophotal shape (on average), but the high diskiness 
of the initial galaxy is washed out by the multiple mi- 
nor mergers. When the total merged mass increases, the 
remnants tend to show boxy isophotes (on average). The 
"orbital scatter" is as in Fig. [5] (see also Sect. I4.1TT|) . 

and sometimes from 2:1 mergers. But these equal-mass 
binary mergers fail to reproduce the most boxy ellip- 
tical galaxies, in particular the giant elliptical galaxies 



tical galaxies, m particular the giant elliptical galax it 
(e.g.. lNaab fc Burkertlliool iNaab fc Ostrikeill200<i ). 



We have measured the 04 parameter for the multiple 
merger remnants from several projected images (same as 
V/a). The average results for the various simulated mass 
ratios are shown in Fig. [5l Once again, the result depends 
mainly on the total merged mass with little dependence 
on how it is merged. For instance, two or three successive 
mergers of mass ratio 7:1 produce the same final low diski- 
ness as one single 3:1 merger, while 10 minor 10:1 mergers 
result in a boxiness comparable, on average, to an equal- 
mass 1:1 merger. Note that, just like a 1:1 remnant, a 
10x10 : 1 remnant can appear slightly disky depending on 
the viewing direction (see the 1-cr scatter in Fig. [5]) but will 
more frequently appear boxy, with on average 04 ^-0.1 - 
-0.2. Nevertheless, a trend can be noted for lower mass 
ratios, in particular 10:1 and 15:1, to converge less rapidly 
towards negative 04 parameters. This indicates that the 
less violent relaxation in minor mergers is less efficient in 
destroying the disky underlying orbital structure. The dif- 
ference remains small after several 10-15:1 mergers, but 
this suggests that very minor mergers at higher mass ratios 
might be less efficient in forming boxy elliptical galaxies. 
This trend could also be increased by the higher gas mass 
fraction in lower mass galaxies, because this gas tends to 
reform disky structures in elliptical galaxies, which can 
reduce the boxiness. 

Our result extends the trend observed in 1:1 vs 3:1 bi- 
nary mergers, namely an increase in the boxiness with the 
merged mass. This trend still holds when several mergers 
are repeated in time. Thus, boxy elliptical galaxies can be 



formed either in the merging of pairs of nearly-equal-mass 
galaxies, but also by the merg i ng of several lower mass 
galaxies. iLima-Neto &: Combed ( 19951 ) had found that re- 
peated mergers tend to wash out any particular shape 
(boxy or disky) of elliptical galaxies. The apparent con- 
tradiction can relate to the much lower resolution available 
to these authors, but also that they studied only collision- 
less and simultaneous major mergers, which is a largely 
different process. 

We actually find here that an increasing number of 
subsequent mergers does increase the boxiness. In the 
same vein, (Naab, Khochfar & Burkert 2006b) find that 
an elliptical galaxy that undergoes another merger (with 
another elliptical galaxy in their simulations) can have 
a slight increase in its boxiness, too. Binary mergers of 
disk galaxies cannot account for the most massive ellipti- 
cal galaxies being boxy: even coUisionless, perfectly equal- 
mass 1:1 mergers generally do not result in very boxy 
systems. The formation of such giant boxy ellipticals can 
then more likely result of from several re-mergers, the first 
ones forming an elliptical galaxy and the last ones mak- 
ing it more and mor e boxy. Th is mechanism is at work 
in the simulations bv lNaab et a l. (2006b), and our results 
(Fig. [5]) further indicate that the sequential mergers do not 
need to be equal-mass nor even major : several 3:1 mergers 
can form a boxy elliptical galaxy, and a sufficient number 
of 7:1 mergers can as well. An alternative scenario to form 
massive boxy elliptical galaxies might be the simultaneous 
merger of several galaxies in a de nse proto-cluster environ- 
ment , as in the simulations by I Weil fc Hernguistl ( 1994 . 
19961 ): the morphological properties of massive elliptical 



galaxies formed this way remain to be studied in detail. 
We at least show that repeated mergers can explain the 
boxiness of large ellipticals, which binary major mergers of 
disks alone fail to explain. It is possible that the giant el- 
lipticals first formed from major mergers at high redshift, 
then continued to grow with repeated minor mergers with 
lower-mass galaxies, hence gradually increasing their box- 
iness. 



3.4. Anisotropy of stellar orbits and the dark matter 
content of elliptical galaxies 

Up to now, wc have analyzed only the projected proper- 
ties of the multiple merger remnants. These can be directly 
measured in observed galaxies. The three-dimensional or- 
bital structure of elliptical galaxies cannot be directly ob- 
served, but there is some evidence that the rotation of 
elliptical galaxi es is generally not large enough to support 
their flattening ( Binnevl[l982 ^. so that their stellar velocity 
dispersion has to be anisotropic. Yet, the actual anisotropy 
of the orbits cannot be accurately inferred from observable 
parameters alone. And indeed, it is only recently that nu- 
merical simulations of major mergers revealed the large 
radial anisotropy of the stellar orbits in elliptical galaxies, 
in particular in the outer regions. This anisotropy results 
in a decrease of the observed velocities (along the line-of- 
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sight) in the outer regions, so that the actual mass content 
of elliptical galaxies had been u nderestimated - in partic- 
ular their dark matter content ( Dekel et al. I I2OO5I) . 



In the multiple merger remnants, we measured the 
anisotropy parameter 



13 = 1- 



where ae and Ur are respectively the tangential and ra- 
dial velocity dispersions with respect to the galaxy mass 
center. The result is shown, on average over all directions 
but as a function of radius, in Fig. [6] The elliptical-like 
galaxies formed in multiple minor mergers tend to have a 
larger radial anisotropy than those formed in one single 
merger bringing the same total mass. 

A possible explanation of the higher radial anisotropy 
found in multiple minor merger remnants can be that 
the stars and gas clouds in smaller companions have a 
lower angular momentum than in a massive one. Then, 
the resulting remnants have a lower angular momentum 
per star, which results in more eccentric orbits, i.e. higher 
radial velocity at a given radius. As a result, higher radial 
anisotropics would te nd to increase th e dynamical mass 
of elliptical galaxies ( Dekel et al. 20051 ) if some of them 
have been formed by multiple minor mergers, but the or- 
der of difference expected from Fig. [6] would not exceed 
a few tens of percent. However, a study including vari- 
ous orbits and initial morphologies is required to confirm 
whether the anisotropy is systemati cally higher in multi- 
ple minor merger remnants. Indeed, lAthanassoula 
shows that many parameters can significantly influence 
the orbital structure of the merger remnant, including the 
time interval between subsequent collisions in group merg- 
ers. Moreover, the radial anisotropy is unlikely to keep 
incre asing with higher nu mbers of minor mergers, given 
that lAthanassoulal (j2005| ) finds lower anisotropy in rem- 
nants whose progenitors are already elliptical. 



4. Discussion 

4.1. Formation of elliptical galaxies without major 
mergers 

4.1.1. Remnants from multiple minor mergers and 
orbital parameters 

Single binary mergers with mass ratios of 5:1 and larger 
do not form remnants with cUiptical-like properties, but 
disk galaxies: these are "minor" mergers, or the "interme- 
diate" mergers forming SO-like galaxies (Bournaud et al. 
2005b). But when several mergers of this type occur, the 
system shows increasing velocity dispersion, together with 
a rounder and more concentrated morphology (lower E 
and higher n). Hence, both the morphological and kine- 
matical properties of these multiple minor merger rem- 
nants are similar to those of observed elliptical galaxies. 
An elliptical-like galaxy is formed when the total mass of 
the companions is 30-40% of the mass of the main initial 




Fig. 6. Anisotropy parameter [3 for an equal-mass major 
merger remnant (1:1) compared to remnants of multiple 
minor mergers (5x5 : 1 and 10x10 : 1). The radial variations 
of the anisotropy are qualitatively similar in these systems, 
but the remnants of multiple minor mergers have a higher 
radial anisotropy in their outer regions. 



galaxy (merged mass = 1.3-1.4). This has been estab- 
lished for 5:1 to 15:1 companions, but the simulation with 
50:1 companions presented below indicates that this result 
extends to smaller companions, too. Multiple minor merg- 
ers are thus a new pathway for the formation of elliptical 
galaxies without major mergers. 

The sequences of mergers in our simulations described 
so far alternate prograde and retrograde companions. The 
impact parameter, velocity, and inclination were not fixed 
at particularly low or high values, and are thus expected 
to give representative results. Furthermore, we know from 
paper I and other studies of binary mergers that the global 
properties of a merger remnant are more influenced by the 
mass ratio of the merging galaxies than the other param- 
eters. For instance, the properties of a binary 3:1 merger 
remnant vary with the orbit, but the variations are gen- 
erally smaller than the typical differences with binary 1:1 
or 5:1 remnants. Yet, the influence of orbital parameters 
may a priori be more important for multiple mergers, be- 
cause the total angular momentum provided to the system 
can change, influencing the residual rotation of the merger 
remnant. 

We have then repeated our 10x10:1 simulation with 
all companions on prograde orbits, all companions on ret- 
rograde orbits, all companions at inclination = (alter- 
nating prograde and retrograde orbits as in the fiducial 
run), and all companions at = 65 degrees. The typ- 
ical variation of the various properties V/a, n and 
tti corresponding to these changes in orbits is shown in 
Figs. [2] to [5l We then notice that: (i) for a given merged 
mass (for instance 10x10:1 vs 3x3:1) the differences be- 
tween the mass ratios are below this orbital scatter, but 
(a) the evolution of the morphological parameters along 
the merger sequences is larger than the scatter related to 
orbital parameters, hence it is a significant general result. 
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Fig. 7. Schematic description of the nature of the binary 
and multiple merger remnants, as a function of the num- 
ber of mergers and the mass ratio. Orbital parameters 
have not been varied in this paper, except for 10:1 merg- 
ers, but we know from paper I that the properties of a 
merger remnant depend on the mass ratio more than the 
orbital parameters, so this scheme should give a repre- 
sentative estimate of the number of mergers required to 
transform a disk galaxy (Sp or SO) into an elliptical-like 
system. The ellipticals formed for instance by 1 : 1 merg- 
ers or 5 X 7:1 sequences are generally boxy on average 
(they appear disky as well under some projections), and 
their boxiness is moderate. An increasing degree of boxi- 
ness can be obtained with an increasing number of merg- 
ers. These results were established for a gas mass fraction 
of 15% within stellar disks but to first order can be applied 
to all systems that arc not initially strongly gas-dominated 
(see text). 



The fact that a 10x10: 1 sequence produces slow- rotating 
and (slightly) boxy ellipticals, on average, while a 6x10 : 1 
sequence results in a disky rotating elliptical is a robust 
result that is not much affected by the variations of orbital 
parameters. Only the case where all the mergers occur on 
prograde orbits gives a somewhat higher V/a, but this is 
a rather unlikely case. In major mergers too, some orbits 
can preserve massive d isk-like structur es (for instance in 
the case of NGC 4550. IPfennige^ll99it . Similarly, a suc- 
cession of ten minor mergers with all companions coplanar 
and corotating to the main galaxy would likely form an 
SO or disky E rather than a boxy E, but such unlikely 
cases can only explain specific situations without being 
representative of the majority of real mergers. 

In Fig. [71 we briefly summarize the properties of the 
binary and multiple merger remnants, as a function of the 
mass ratio and number of mergers, according to the sim- 
ulations of this paper and that of paper I. Real situations 
can obviously be more complex, but from what we said 
above, we can expect for instance that a 5 : 1-10 : 1-10 : 1 
sequence will form an clliptical-like galaxy whose average 
properties are similar to the 4x10: 1 and 2x5: 1 remnants. 



4.1.2. Higher mass ratios: can very small companions 
still form ellipticals? 

Here we present a simulation where the same initial spi- 
ral galaxy as before merges with 50 companions each of 
l/50th of its initial mass. The parameters are the same as 
for the other simulations, except that: 

- the mergers cannot be fully subsequent; they are uni- 
formly distributed over a 8 Gyr period. 

- the dwarf companions initially contain 30% of their vis- 
ible mass in gas. 

The relaxed remnant of this 50x50 : 1 multiple minor 
merger sequence is displayed in Fig. [8l viewed along the 
line-of-sight giving the largest apparent flattening (which 
is E = 6.1). This system has the shape of an E6 elliptical 
galaxy (but would more likely be observed as an E5 to EO 
one under other projections), and is dominated by veloc- 
ity dispersions. We show the mean velocity and velocity 
dispersion profiles that would be observed along the ma- 
jor axis of this system in Fig. [21 which confirms that it 
has elliptical-like large-scale kinematical properties with 
V/a 0.5 even in the outer regions. 

The remnants of 5:1 to 15:1 multiple mergers did not 
show a large variation of their properties with the mass ra- 
tio at given total merged mass. This already suggests that 
the result was the same for a higher number of smaller 
companions. The result of the 50x50: 1 confirms that the 
formation of elliptical-like galaxies can still be achieved 
with high mass ratios. In detail the 50x50 : 1 merger rem- 
nant is slightly fiatter and has a somewhat higher residual 
rotation than the other cases; the discrepancy is not major 
and likely results from the higher gas fraction in the 50:1 
companions. Thus, a large number of mergers with very 
high mass ratios still form elliptical- like galaxies, at least 
when the size/concentration of galaxies is scaled as in our 
model (which assumes M oc E? and a constant central 
surface density). Resolution tests should be performed in 
the future to further explore the case of such very small 
companions. A few 10^ particles per galaxy is usually be- 
lieved to provide a viable large-scale description for major 
mergers, but in the 50:1 cases the small companion galaxy 
contains only 2 x 10^ stellar particles. High mass ratios are 
difficult to model for this reason, and a lack of resolution 
could for instance smooth their gravitational potential and 
underestimate their effects. 



4.1.3. Mergers of gas-rich and high-redshift galaxies 

The initial gas mass fraction in the initial spiral galaxy 
and in the 1:1 to 15:1 companions is 15% in our models: 
the mergers are not "dry" , but higher gas fractions could 
still be encountered in real mergers at high redshift and/or 
in low-mass companions. A high fraction of gas in mergers 
tends to preserve the diskiness and somewhat higher resid- 
ual rotation (Naab, Jesseit & Burkert 2006a; Jesseit et al. 
2007) because the gas can be stripped before the galaxy 
collision, making it less violent, and also because gas can 
fall back into a disk where it will form new stars after 
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Fig. 8. Projected mass density, veloc- 
ity and dispersion fields for a relaxed 
50x50:1 merger remnant, viewed along 
the projection resulting in the largest 
apparent flattening {E — 6.1). The con- 
tours on the velocity field are 35 and 70 
km s~^, on the velocity dispersion field 
they are 100, 120 and 140 km s"^ Note 
the misalignment between the residual 
rotation and apparent flattening axis. 
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Fig. 9. Average velocity and velocity dispersion for the 
50x50 : 1 merger remnant shown in Fig. [HI taken along a 1 
kpc-width slit aligned with the apparent major axis. The 
system is supported by large velocity dispersions similar 
to those observed in elliptical galaxies. 



ing more than half of the mass, will to first order end-up 
in the spheroidal component after a sufficient number of 
minor mergers. Only the stars born during and after the 
mergers potentially end up in a disk, but not all of them 
(for instance not those formed during the first merger in 
a sequence of ten mergers), so that the disk component 
cannot dominate the mass distribution - which it does in 
spiral galaxies. Thus, at redshift z ~ 1 and below, we can 
predict that a succession of minor mergers with the mass 
ratios studied so far will still transform spiral galaxies into 
ellipticals, even if the detailed properties of the latter can 
vary somewhat with the gas fractions, and the exact num- 
ber of mergers required to make an elliptical could change 
too. Note though that high-redshift disks frequently have 
peculiar morphology with sometimes much mass in giant 
clumps (Conselice 2003; Elmegreen, Elmegreen & Hirst 
2004; Conselice et al. 2004; Elmegreen 2007; Bournaud, 
Elmegreen & Elmegreen 2007), and the potential infiu- 
ence of such properties on the structure of high-redshift 
major and minor merger remnants remains to be explored. 



the merger. However, major mergers rarely preserve disk 
galaxies similar to real spirals. Indeed, even cases start- 
ing with pure gas disks in a coplanar co-rotating situation 
(the best case to preserve most of the angular momentum) 
end up with a bout half the mass in a sphe roidal compo- 
nent (see e.g., 'Springel fc Hernguist 2005f) . More realis- 
tic orbits will have more mass in the spheroid, and more 
realistic initial gas fractions will reduce the subsequent 
disk re-formation so that even more mass will end up in 
a spheroid, making the galaxy elliptical-like (there can 
be an outer disk component, formed for instance by gas 
falling-back, but not dominating the mass). So gas- rich 
major mergers are still expected to form elliptical galax- 
ies, even if their detailed properties show differences. Only 
extreme situations with very high initial gas fractions can 



preser ve an (early-type) disk (see also ,Robertson et al. 



( 2006^ ). which can be representative only of primordial 
galaxies at very high redshift, but not of most galaxies at 
z ^ 1. Our simulations, even with 30% gas for the 50:1 
case, can have a low-mass disk component but not a final 
mass distribution dominated by an exponential disk. 

As for multiple minor mergers, except at very high red- 
shift z >> 1, gas mass fractions are not expected to be 
larger than 50%, except perhaps in some exceptional cases 
or in LSB galaxies. Then the pre-existing stars, represent- 



4.1.4. Comparison with major mergers and star 
formation history 

The global properties of the multiple minor merger rem- 
nants resemble those of binary major mergers. It is thus 
difficult to infer from the projected properties of an ellip- 
tical galaxy its past formation mechanism. The differences 
between multiple minor and major mergers are not larger 
than the scatter related to orbits and other parameters. 
The orbital structure can be different, but it is unclear 
whether there is a systematic difference, because many 
parameters that may also influence the anisotropy of the 
merger remnants have not been varied in this paper. 

Differences may more likely be found using more de- 
tailed parameters, like the /13 of the projected velocity dis- 
tribution (Gonzalez-Garcia, Balcells fc Olshevsky 2006), 
the Xfi parameter defined bv lEmsellem et all (|2007t ) which 
couples the resolved kinematics to the density profile and 
revealed different families of early- type galaxies, or corre- 
lations of several such parameters. Such detailed compar- 
isons with advanced parameters are beyond the scope of 
the present paper which studies whether merger remnants 
are elliptical-like or not, but will be important for future 
studies. 
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Fig. 10. Star formation history in binary major mergers 
(1:1, for the prograde and retrograde orbits) and for a 
succession of 7:1 minor mergers. 
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Another major difference lies in the star formation 
history of these systems. Ind eed, minor mergers trigger 
bursts of star formation (e.g.. ICoxll2004t ICox et al.ll2007( ) 
that ar e comparative ly less intense than in major merg- 
ers (see lDi Matteo et al. 2007a, regarding the intensity of 
starbursts in major mergers). As a result, the star for- 
mation history of multiple minor merger remnants largely 
differs from that of a binary major merger, as illustrated 
in one case in Fig. 1101 This difference may however be- 
come smaller if the successive minor mergers occurred 
with smaller time intervals or simultaneously. 

4.1.5. Scaling relations in elliptical galaxies 

Beyond their fundamental properties studied above (mass 
profile, kinematics, isophotal shape), the real ellipti- 
cal galaxies are also observed to be described by rela- 
tions between these fundamental parameters. An impor- 
tant sca ling relation is the so-called fundamental plane 
(jFaber et al., 1987 ). Most scaling relations are defined over 
the complete mass (or luminosity) range of elliptical galax- 
ies, while our simulations explore only a restricted mass 
range (a factor of two between the initial and final masses), 
making it irrelevant to test the viability of multiple minor 
mergers by a comparison with these relations. This is also 
the case for major merger simulations, and simulations 
covering the full mass range of ellipticals in a cosmologi- 
cal context would be required. 

An observed fundamental relation of ellipticals that 
can be tested in our simulation s is the ellipti city - veloc- 
ity dispersion [E, V/cr) relation ( Binnevl[l982 h. which does 
not explicitly span a large mass range. In our simulations 
of multiple minor mergers, we find in Fig.[TT]that a rather 
tight relation exists between E and V/a for the elliptical 
remnants of multiple minor mergers - and even for those 
that we classify as disks (spirals or SOs) after for instance 
one single 7:1 merger or two 15:1 mergers. The multiple 
merger mechanism then appears to form elliptical galax- 



Fig. 11. Relation between the velocity dispersion (mea- 
sured via the V/a parameter) and ellipticity (flattening E) 
along the multiple minor merger sequences with various 
mass ratios, and the reference major mergers. Evolution 
with increasing number of mergers is from top-right to 
bottom- left. 



ies that are viable with regard to this observed relation. 
The fact that all fundamental parameters have a similar 
evolution along all mergers sequences, whatever the mass 
ratio is (Figs. 2 to 5), suggests that other parametric re- 
lations exist too. Note however that we cannot check in 
our simulations whether higher- or lower-mass ellipticals 
would still lie on the relation shown in Fig. [TT] and not 
offset from it. 

The mass of spheroids (bulges and ellipticals), of- 
ten traced by central velocity dispersion, is known to 
correlate with the mass of central black holes (e.g., 
Kormendv fc Richstonel Il995l ). probably mostly grown 
durin g rapid accretion QSO phases (jYu fc Tremaine 
20021 ). Major mergers c an directly fuel a centra, l black hole 
by rapid gas accretion (jPi Matteo et al.ll2007br ). while mi- 
nor mergers would rather lead to more moderate inflows 
(see the star formation history in major vs multiple mi- 
nor mergers in Fig. I10|) . Still, minor mergers can increase 
the mass of a central black hole by the merging of several 
smaller black holes - those initially at the center of the 
primary target galaxy and each companion: then the rela- 
tion between the central black hole mass and the spheroid 
mass could a priori be preserved in the elliptical galaxies 
formed by multiple minor mergers. 



4.2. Constraints on galaxy growth and evolution 

Repeated minor mergers, even with rather high mass ra- 
tios, are thus a mechanism that can form spheroids hav- 
ing the properties of the real observed elliptical galaxies. 
Multiple minor mergers is then a possible mechanism to 
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form elliptical galaxies without major mergers, at least 
within the theoretical frame where binary and repeated 
mergers plus internal evolution are the only phenomenon 
driving the evolution of galaxies. The real evolution of 
galaxies is more complex, because galaxy-galaxy mergers 
and internal evolution are not the only driving mecha- 
nisms, and we now discuss how repeated minor mergers 
can be placed in a more complete description of galaxy 
evolution. 



4.2.1. The survival and growth of spiral disks 

Disk galaxies are observed to increase their size and 
mass between z = 1 and z — 0. At redshift 1, 
stellar disks had on average smaller truncation radii 



(jTrujillo &: Pohlenll2005r) and smaller exponential scale- 
lengths (e.g., lElmegreen et al. l2007bl) . Their mass must 
then have increased in significant proportions over the last 
~ 8 Gyr. If this was to be achieved mainly through minor 
mergers, most spiral galaxies would have been transformed 
into SO (mass increased by ~20-30% through mergers) 
and ellipticals (mass increased by 30-40% and more). It 
is then required that another process participates in the 
growth of spiral galaxies, likely the accretion of external 
gas. While minor mergers progressively destroy disks, ac- 
cretion of gas can refuel massive thin disks to maintain 
high disk-to-bulge mass ratios. It can be the accretion of 
diffuse gas falling into the potential well of galaxy ha- 
los, or the accretion of gas-dominated clumps with low 
visible and dark masses, that are disrupted in the early 
stages of the interaction. But usual dwarf galaxies of 
Mvis ~ 10^ Mq and above plus their massive dark mat- 
ter halo cannot achieve this, as shown for instance by our 
simulation with 50:1 mergers. 

Cold gas can be accreted directly from cosmological fil- 
aments for galaxies whose total mass (including the dark 
halo) does not exceed ^ 10 ^^ Mq IjBirnboim fc Dekell 
2003tlDekel fc Birnboimll2006f) and even for more massive 



galaxies but less regularly (Birnboim, Dekel & Neistein 
2007); and it can explain the morp hological properties of 
disk gal axies in our Local Uni verse ( Bournaud fc Combed 
[2002; Bournaud et al. 2005a ). Moreov er, the growth of 
disks is observed to be inside-out (e.g., Truiillo fc PohlenI 
2005h . which cannot be easily explained by a (minor) 



merger-driven growth, but can result from the accretion 
of gas if the lower angular momentum material is accreted 
first, and higher angular momentum gas accreted at later 
stages. Gas falling-back from large radii after the merg- 
ers (in particular gas from tidal tails) appears not to be 
massive enough to refuel a massive disk, even in our 50:1 
cases with gas-rich companions. A low-mass disk compo- 
nent can be fueled this way, but a disk dominating the 
total mass distribution cannot be maintained without the 
additional accretion of external gas. 

That spiral disks have been able to grow from 2; = 1 to 
z = thus implies that some phenomenon favor the per- 
sistence of their disks, like the accretion of cold gas. Then, 



the process of progressive disk destruction and transfor- 
mation into ellipticals by repeated minor mergers could be 
compensated for by other processes, like cold gas accre- 
tion, which would increase or maintain massive disk com- 
ponents around/within the heated stellar spheroid. But in 
our simulations the galaxy evolution is driven only by mi- 
nor mergers, and this leads to a progressive transformation 
of disk galaxies into elliptical- like spheroids. 

4.2.2. The number of elliptical galaxies and their 
formation 

The standard "major merger" scenario for the formation 
of elliptical galaxies is usually assumed to produce a frac- 
tion of elliptical galaxies roughly consistent with that ob- 
served (e.g., Mamon 1992; Baugh, Cole fc Fren k 1996). An 
additi onal mechanism has been proposed by Naab et al.l 
(|2007l ). but is mainly expected to be at work in the Early 
Universe forming primordial elliptical galaxies. 

Successions of several minor mergers are mor e likel y 



see 



Khochfar fc sil^ (|2006l) . 



than binary major mergers , 

iMaller et al. ( 2006h . and estimates in paper I - at least 
at redshift z < 1. The number of elliptical galaxies is 
then a priori expected to increase significantly through 
this alternative mechanism. But then the number of el- 
liptical galaxies present in the Universe today would be 
much larger than observed. Even among field spiral galax- 
ies, which have significantly grown over the past 8 billion 
years, many would have become elliptical and lenticular 
galaxies. Another process is then required to prevent the 
succession of minor mergers from systematically trans- 
forming disk galaxies into ellipticals as soon as their mass 
is increased by a few ten of percents, otherwise this mech- 
anism together with binary major mergers would result 
in an excess of elliptical galaxies. Still, there should be at 
least some elli pticals formed by m ultiple minor mergers 
at low redshift; van Dokkum ( 20051 ) finds a large fraction 
of mass ratios larger than 4:1 in the mergers undergone 
recently by field ellipticals. 

Observationally, on-going minor mergers are not un- 
common at z = 0, but also at z ~ 1: for instance a large 
fraction of the interacting systems at z ^ 1 in GE MS and 
GOODS is of M51-type (|Elmegreen et al.ll2007a[ l. which 
corresponds to mergers where the mass of the companion 
is far from negligible, as in our 15:1 - 5:1 mergers. One 
merger of this type can transform a spiral into an Sa or 
SO (paper I) and a few ones make it elliptical. Galaxies 
thus have to accrete cold gas in a diffuse phase, or small 
clumps that are not too massive/concentrated so that they 
are disrupted before they reach the disk. 

Nevertheless, multiple minor mergers must have 
formed some elliptical galaxies. In particular, some 
galax ies that are too massive to accrete cold gas 
(Bi rnboim fc Dekell |2003[ ) or in dense environments have 
likely evolved into ellipticals under the combined ef- 
fect of repeated minor mergers (in addition to binary 
major mergers). This provides a good explanation for 
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the boxiness of giant ellipticals, as discussed previously 
(Section 3.3), and there is further evidence that ellipti- 
cal galaxies did not all fo rm from binary major merg- 
ers (jNaab fc Ostrikerj|2007l) . Repeated major mergers in 
groups may have formed them too, but these are less likely 
than multiple minor mergers. Thus, successions of minor 
mergers likely participated in the formation of some of the 
elliptical galaxies seen today. 

Elliptical galaxies have a higher frequency of globular 
clusters ( GCs) relatiye to their stellar mas s than spiral 
galaxies ()Harrisl 1991 ; van den Bergh 2001 ) . The forma- 
tion of GCs can be triggered b y strong shocks associated 
with starbu rsts, as suggested bvlvan den Berghl(|l979f ) and 
studied by |Jog fc Solomon ( 1992 ). Major mergers trig- 
ger stronger starbursts, while the shocks and starbursts 
are less intense (but repeated) during the successive mi- 
nor mergers (see Sect. 14.1^ . potentially forming different 
numbers of GCs - which models resolving GC formation 
could confirm. This could be a way to disentangle the 
elliptical galaxies formed by major and repeated minor 
mergers. 

5. Conclusion 

In this paper, we have shown that the succession of minor 
galaxy mergers (mass ratios larger than 4:1) leads to the 
gradual transformation of spiral galaxies into elliptical-like 
galaxies. This is the case both for the truly minor mergers 
(mass ratios above 10:1, a single merger of this type forms 
a spiral) and the "intermediate" mergers (mass ratios be- 
tween 4:1 and 10:1, a single merger of this type forms an 
SO- like galaxy, e.g. paper I). The remnants of repeated 
mergers with these mass ratios have a spheroidal shape, a 
radial profile close to the "i^i/^" empirical profile of ellipti- 
cals, and are supported by large velocity dispersions V ja. 
This is true even for very low-mass companions (up to 
50:1 in one of our simulations), and the relaxed remnant 
resembles an elliptical galaxy as soon as the total mass 
added by successive mergers exceeds 30-40% of the mass 
of the main initial spiral. The global properties of the mul- 
tiple merger remnants depend more on the total merged 
mass than on the mass ratio of each merger, and for in- 
stance a 10:10x1 merger sequence resembles a 1:1 rem- 
nant more than a 5 : 10x1. These properties can vary with 
the orbital parameters and morphology of the progenitor 
galaxies, but not widely. In particular, only extremely gas- 
rich systems, or companions on peculiar orbits, may keep 
disk-dominated galaxies after several minor mergers, but 
this should not be the most frequent case except at very 
high redshift. A noteworthy dynamical result is that the 
final properties of the remnant (radial profile, flattening, 
kinematics) are to first order independent of whether it 
is formed by a single 1 : 1 merger, a 4 x 4:1 sequence, or 
a 10 x 10 : 1 sequence, so as long as the final mass has 
increased by the same factor. 

Multiple, sequential minor mergers provide a new the- 
oretical pathway to form elliptical galaxies without major 
mergers. In a purely hierarchical scenario, a succession of 



several minor mergers is more likely than one single binary 
merger, and must then have formed some of the elliptical 
galaxies seen today. Repeated minor mergers also provide 
a possible explanation for some specific features, in par- 
ticular the boxiness of the largest ellipticals, that are not 
explained by the standard major mergers (1:1-3:1). 

However, minor mergers repeated in time would tend 
by themselves to produce too many elliptical galaxies. 
Moreover if field galaxies were to grow only by hierarchical 
merging, most spirals would have been transformed into 
ellipticals over the last 8 Gyr, even if these mergers were 
minor ones. It is then required that another process par- 
ticipates in the growth of galaxies, braking down or revers- 
ing the minor-merger-driven evolution towards early-type 
spheroids. This process is likely to be accretion of cold gas, 
either in a completely diffuse phase, or in small gas-rich 
clumps - but not dwarf galaxies bounded by their dark 
matter haloes that would have destroyed massive disks. 

While we present a new mechanism for the formation 
of elliptical galaxies in a purely hierarchical context, we 
also conclude that its efficiency must be somewhat limited: 
the evolution of galaxies cannot be purely merger-driven, 
and a large fraction of their mass cannot have been ac- 
quired through mergers, even minor ones with small com- 
panions. The growth of spiral galaxies at low and moder- 
ate redshifts must have rather been a competition between 
these minor mergers and accretion of cold gas. Some of 
the present-day ellipticals likely formed by multiple minor 
mergers, but this process cannot be much more efficient 
than binary major mergers, otherwise it would result in 
an excess of elliptical galaxies. 

Cosmological simulations of restricted volumes can 
reach a resolution similar to that of the generic galaxy 
simulations pres ented in this paper (see in particular 
Naab et al.ll2007l ). but not over large samples that can be 
used for statistical purposes. At the opposite, large- volume 
cosmological simulations may lack the resolution needed 
to reproduce the effects of repeated minor mergers in a re- 
alistic way, because this requires that the dwarf compan- 
ions are sufficiently resolved. However, these simulations 
can provide accurate predictions of the evolution history 
of galaxies (mass ratios of mergers, time intervals, and 
rate of diffuse gas accretion). The possibility of forming 
ellipticals but keeping most field galaxies as spirals at the 
same time, in high-resolution galaxy models taking into 
account these large-scale cosmological predictions, could 
be a future test of cosmological models. 
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